Comment on "Numerical estimates of the spectrum for anharmonic PT symmetric 

potentials" by Bowen et al 
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The paper by Bowen, Mancini, Fessatidis, and Murawski (2012 Phys. Scr. 85, 065005) demon- 
strates in a dramatic fashion the serious difficulties that can arise when one rushes to perform 
numerical studies before understanding the physics and mathematics of the problem at hand and 
without understanding the limitations of the numerical methods used. Based on their flawed nu- 
merical work, the authors conclude that the work of Bender and Boettcher is wrong even though it 
has been verified at a completely rigorous level. Unfortunately, the numerical procedures performed 
and described in the paper by Bowen et al are incorrectly applied and wrongly interpreted. 
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The key sentence in the paper by Bowen et al [l[ is in 
Section 3 (Discussion) : "The mystery of this study is why 
does this numerical method give such radically different 
results from those discussed by Bender." As we explain 
in the present Comment, the simple and straightforward 
explanation of this "mystery" is that the results reported 
in Bender and Boettcher [2j are correct and that the nu- 
merical work reported in the paper by Bowen et al is 
incorrect and misinterpreted. 

In the original paper by Bender and Boettcher Q , the 
eigenvalues of the "PT-symmetric Hamiltonian 



H = p 2 + x 2 {ixf, 



(1) 



where e is a real parameter, were studied in great detail 
by using both numerical and perturbative methods. It 
was shown that when e > 0, the eigenvalues are real, pos- 
itive, and discrete and that there are no complex eigen- 
values. These conclusions were subsequently verified in 
rigorous studies by Dorey, Dunning, and Tateo [3j and 
Shin % 

Before one can perform numerical calculations of the 
eigenvalues of the Hamiltonian (JT|), it is essential to un- 
derstand what it means to extend an eigenvalue problem 
into the complex domain. The first detailed work in this 
area was done by Bender and Wu Q , who examined the 
behavior of the eigenvalues of the quartic anharmonic os- 
cillator as functions of complex coupling constant. Fur- 
ther studies of analytically continued eigenvalue prob- 
lems were done by Bender and Turbiner |6|. The crucial 
element explained in all of these studies is that as a pa- 
rameter in an eigenvalue differential equation is varied, 
the Stokes wedges in which the boundary conditions on 
the eigenfunctions are imposed must rotate as functions 
of that parameter. 

To be specific, the eigenvalues associated with the 
Hamiltonian (JT|) are determined by the complex differ- 
ential equation 



ip"(x) + x 2 (ix) e %k{x) = Eif>(x) 



(2) 



and a pair of PT-symmetric boundary conditions in the 
complex- x plane: ip(x) — > as \x\ — > oo with argx lying 
inside the opening angles of the Stokes wedges. 

The details of how to find the opening angles of Stokes 
wedges are not given here because detailed explanations 
are given in Refs. 0, 0, 0- We simply state the results: 
When e = 0, the two Stokes wedges have opening an- 
gles of 90° and are centered about the positive-real and 
negative-real axes. As e increases from 0, the Stokes 
wedges become thinner and rotate downward into the 
lower-half complex- a; plane. When e reaches the value 
2, the opening angles of the wedges have decreased to 
60° and the wedges no longer contain the positive- and 
negative-real axes. 

Because the eigenvalue differential equation @ must 
be solved along a contour in the complex- x plane that 
terminates in the Stokes wedges, it is best to solve the 
differential equation numerically on a complex contour 
by using the standard Runge-Kutta method. The Runge- 
Kutta procedure described in Refs. [U, H| has the obvious 
advantage that it can be used for any real value of e, 
and not just for integer e. It is also extremely important 
to understand that different pairs of Stokes wedges give 
different sets of eigenvalues. For example, for the quan- 
tum harmonic oscillator Hamiltonian H — p 2 + x 2 , which 
corresponds to e = in ([1]), if Stokes wedges containing 
the real axes are used, the eigenvalues are 1, 3, 5, 7, . . ., 
which are strictly positive. However, if Stokes wedges 
containing the imaginary axes are used, the eigenvalues 
are —1, —3, —5, —7, . . ., which are strictly negative. 

In the paper by Bowen et al a clumsy and inappropri- 
ate numerical procedure for calculating the eigenvalues 
is used in which the Hamiltonian is expanded in terms 
of Harmonic-oscillator basis functions. The procedure is 
clumsy because the algebra becomes much too unwieldy 
when e is noninteger. Thus, in the paper by Bowen et 
al the procedure was limited to the very special cases 
e = 0, 1, 2, 4, 6. 

Furthermore, while the harmonic-oscillator eigenfunc- 
tions are complete on the real axis, these basis functions 



2 



are not complete in the complex plane, so these basis 
functions cannot be used unless the Stokes wedges in- 
clude the real axis.^ Thus, the numerical calculations 
that were performed by Bowen et al are completely in- 
valid and meaningless except for the trivial case of the 
harmonic-oscillator (e — 0) and the more interesting and 
nontrivial case e = 1. 

Unfortunately, Bowen et al made a long sequence of 
wrong arguments and completely misinterpreted the nu- 
merical results that they obtained for the e = 1 case. 
While it is perfectly correct to expand H = p 2 + ix 3 in 
an infinite harmonic-oscillator basis, Bowen et al then 
truncate the infinite matrix representation of the Hamil- 
tonian, which is VT symmetric, to an N x N matrix, 
which is no longer VT symmetric. As a result, the finite- 
dimensional matrix has complex eigenvalues. These com- 
plex eigenvalues reported by Bowen et al are merely ar- 
tifacts of the truncation procedure that they used. 

Bowen et al then go on to make further serious mis- 
judgments. The procedure of truncating an infinite ma- 
trix is variational in character. Thus, it can only be used 
to compute the low-lying eigenvalues. If the procedure 
is used carefully and properly, one follows the behav- 
ior of the low-lying eigenvalues as N increases and ob- 
serves that at first the behavior is irregular and complex. 
(The behavior of the eigenvalues is not monotone because 
the N x N matrix is not Hermitian.) However, as N 
gets larger, one observes that the eigenvalues settle down 



one-by-one starting at the low-energy end of the spec- 
trum, and stabilize at the real values found by Bender 
and Boettcher. This stabilization process is extremely 
slow; for very large N, say 100, only about half a dozen 
eigenvalues can be determined with any useful accuracy. 
The rest of the eigenvalues of the N x N matrix bear no 
resemblance to the true eigenvalues of the Hamiltonian 
H = p + ix 3 . Unfortunately, Bowen et al take their nu- 
merical results for the high-lying eigenvalues seriously in 
their paper even though they have no numerical accuracy 
whatsoever. Indeed for the p 2 + ix 3 Hamiltonian, WKB 
predicts that the nth eigenvalue grows like n 6 ^ 5 for large 
n, whereas the numerical work of Bowen al gives nothing 
of the sort. 

Before concluding, we point out that Bowen et al could 
have benefited by actually reading and understanding 
the review paper Ref. Q- They clearly did not do so 
because no less than five times in their six-page paper 
they make the ridiculous claim that in this review pa- 
per Bender states that the Hamiltonian H = p 2 ~~ x 2 has 
real, discrete, negative eigenvalues. In fact, this Hamil- 
tonian certainly does not have real negative eigenvalues, 
and Bender has never made such an absurd claim in any 
paper. 
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